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ABSTRACT 

In order to examine how narrow emission-line flux ratios depend on the Seyfert type, we compiled 
various narrow emission-line flux ratios of 355 Seyfert galaxies from the literature. We present in this 
paper that the intensity of the high-ionization emission lines, [Fe vn]A6087, [Fe x]A6374 and [Ne v]A3426, 
tend to be stronger in Seyfert 1 galaxies than in Seyfert 2 galaxies. In addition to these lines, [0 iii]A4363 
and [Ne iii]A3869, whose ionization potentials are not high (< 100 eV), but whose critical densities are 
significantly high (;> 10 7 cm -3 ), also exhibit the same tendency. On the other hand, the emission-line 
flux ratios among low-ionization emission lines do not show such a tendency. We point out that the 
most plausible interpretation of these results is that the high-ionization emission lines arise mainly from 
highly-ionized, dense gas clouds, which are located very close to nuclei, and thus can be hidden by 
dusty tori. To examine the physical properties of these highly-ionized dense gas clouds, photoionization 
model calculations were performed. As a result, we find that the hydrogen density and the ionization 
parameter of these highly- ionized dense gas clouds are constrained to be nn > 10 6 cm -3 and U > 10 -2 , 
respectively. These lower limits are almost independent both from the metallicity of gas clouds and from 
the spectral energy distribution of the nuclear ionizing radiation. 

Subject headings: galaxies: active - galaxies: nuclei - galaxies: quasars: emission lines - galaxies: 
quasars: general - galaxies: Seyfert 



1. INTRODUCTION 

Seyfert nuclei are typical active galactic nuclei (AGNs) 
in the nearby universe. They have been broadly classified 
into two types based on the presence or absence of broad 
(typically ;> 2000 km s -1 ) permitted emission lines in their 
optical spectra (Khachikian, Weedman 1974); Seyfert nu- 
clei with broad lines are type 1 (hereafter SI), while those 
without broad lines are type 2 (S2). This difference is 
thought to be due to a dependence of the visibility of 
broad-line regions (BLRs) on the viewing angle. Accord- 
ingly, the BLR is thought to be located in a very inner re- 
gion (a typical radial distance from the central black hole 
is r ~ 0.01 pc; see, e.g., Peterson 1993) and surrounded 
by a geometrically- and optically-thick dusty torus (AGN 
unified model; see Antonucci 1993 for a review). 

Contrary to the BLR emission, narrow (typically < 1000 
km s _1 ) permitted and forbidden emission lines, arising 
from narrow- line regions (NLRs), are seen in the spec- 
tra of both Sis and S2s. Therefore, the NLR is believed 
to be located far from the nucleus, and thus not to be 
hidden by dusty tori. If this is the case, the observed 
physical properties of ionized gas in NLRs do not depend 
on the viewing angle toward dusty tori. However, several 
studies have statistically shown that some high-ionization 
emission lines, such as [Fe vn]A6087, [Fe x]A6374, and [Ne 
v]A3426, are stronger in spectra of Sis than in those of S2s 
(e.g., Shuder, Osterbrock 1981; Cohen 1983; Murayama, 
Taniguchi 1998a; Schmitt 1998; Nagao et al. 2000). Some 
emission lines whose ionization potential is not so high 
(< 100 eV), but whose critical density is high, such as [O 
iii]A4363 and [Ne iii]A3869, also show the same Seyfert- 



type dependence (e.g., Osterbrock et al. 1976; Heckman, 
Balick 1979; Shuder, Osterbrock 1981; Schmitt 1998; Na- 
gao et al. 2001b). These Seyfert-type dependences in some 
emission-line intensities seem to conflict with the frame- 
work of the AGN unified model. Therefore, some possible 
models have been proposed to explain such Seyfert-type 
dependences in the narrow emission-line intensities. One 
of the proposed ideas is that these Seyfert-type depen- 
dences in the strength of the high-ionization emission lines 
are caused by a viewing-angle dependence of the visibil- 
ity of highly-ionized dense gas clouds, which are located 
very close to the nucleus, and thus can be hidden by dusty 
tori (Murayama, Taniguchi 1998a, b; Nagao et al. 2000, 
2001b). On the contrary, Schmitt (1998) pointed out that 
the Seyfert-type dependences in the intensities of the high- 
critical-density transition can be understood if the intrinsic 
difference in the NLR size between Sis and S2s (see, e.g., 
Schmitt, Kinney 1996) is taken into account. 

Which is the origin of the Seyfert-type dependence of 
those emission lines, the inclination effect (i.e., obscura- 
tion by dusty tori) or the intrinsic difference in the size 
of NLRs? To investigate this issue, it should be examined 
how emission-line strengths depend on the Seyfert-type. 
Thus, we have compiled emission-line flux ratios of many 
Seyfert galaxies from the literature. In this paper, we show 
the Seyfert-type dependence of various emission-line flux 
ratios based on this compiled database and discuss the 
origin of the Seyfert-type dependence. 

2. THE DATA 
2.1. Data 
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In this section, we briefly describe the properties of the 
database of emission-line flux ratios. See Nagao (2001) 
for details of this database. The database contains vari- 
ous emission-line flux ratios of 355 Seyfert nuclei in total; 
33 narrow-line Sis (NLSls; see, e.g., Osterbrock, Pogge 
1985), 48 broad-line Sis (BLSls), 48 Seyfert 1.2 galax- 
ies (S1.2s), 78 Seyfert 1.5 galaxies (SI. 5s), 9 Seyfert 1.8 
galaxies (SI. 8s), 28 Seyfert 1.9 galaxies (S1.9s), 6 S2s with 
broad emission in near- infrared spectra (S2nir-blrs), 12 
S2s with broad polarized emission (S2hblrs), and 93 S2s 
without any broad emission in their spectra (S2 _ s). In this 
paper, SI. 2s are basically included in BLSls, and SI. 8s, 
SI. 9s, and S2nir-blrs are gathered into the class of "S2 
with reddened BLR (S2rblr)"- When necessary, NLSls 
and BLSls are gathered into the class of "Sl to tai", and 
S2rblr, S2hblrs, and S2~s are gathered into the class of 
"S2 to tai" ■ We adopt the classification of the Seyfert types 
by Veron-Cetty and Veron (2000). 

Since it is often difficult to measure the flux of narrow 
Balmcr components accurately for Sis (and SI. 5s), the 
correction for dust extinction adopting the Balmer decre- 
ment method (see, e.g., Osterbrock 1989) would cause pos- 
sible systematic errors. Therefore, we did not correct any 
dust extinction effect on the emission-line flux ratios in the 
database. The effects of dust extinction on the following 
discussion are mentioned in subsection 3.2. 

2.2. Selection Biases 

The sample used here is not a complete one in any sense. 
Therefore, it is necessary to check whether or not the sam- 
ple is appropriate for the statistical and comparative in- 
vestigations carried out in the following sections. Since 
some possible biases may be caused if there are large sys- 
tematic differences in their redshift or intrinsic luminosity 
distributions, we check their distributions below. 

First, we check the frequency distribution of redshift for 
each type of Seyfert galaxies. The histogram of the fre- 
quency distribution, the median, the average, and the 1 a 
deviation of the redshift for each type of Seyfert galax- 
ies arc presented in figure 1 and table 1. It appears 
that the Sis and the SI. 5s tend to have higher redshifts 
than the S2s. In order to check whether or not the fre- 
quency distributions of the redshift are statistically differ- 
ent among the various types of Seyfert galaxies, we apply 
the Kolmogorov-Smirnov (KS) statistical test (see, e.g., 
Press et al. 1988). The null hypothesis is that the red- 
shift distributions of two types of Seyfert galaxies (in any 
combination) come from the same underlying population. 
The derived KS probabilities (i.e., the probabilities of two 
samples drawn from the same parent population) are sum- 
marized in table 2. The KS test leads to the following 
results: (1) the redshift distributions of the NLSls, the 
BLSls, and the SI. 5s are statistically indistinguishable, 
(2) those of the S2rblr, S2hblr, and S2~ are also statis- 
tically indistinguishable, while (3) the NLSls, the BLSls, 
and the SI. 5s have systematically higher redshifts than 
the S2rblr, S2hblr, and S2~. Therefore, we must keep 
in mind that some properties investigated in the follow- 
ing sections may be affected by the redshift bias. We will 
check this possibility when we carry out statistical investi- 
gations of the emission-line intensity ratios (see subsection 

1 In this paper, a cosmology of Hg = 50 km s~ x Mpc -1 and qo = 
Moshir et al. (1992). 



3.2). 

Second, we check the intrinsic luminosity distribution of 
each type of Seyfert galaxies. Following the AGN unified 
model, the nuclear nonthermal continuum radiation of S2s 
is absorbed by dusty tori and cannot be observed directly. 
This results in weaker continuum emission in the sample 
of S2s than in the sample of Sis, if the distribution of 
the intrinsic luminosity is the same between the samples 
of Sis and S2s. In other words, we may pick up intrinsi- 
cally luminous S2s compared to Sis by survey observations 
with a certain limiting flux. Since such biases may affect 
the statistical properties of our samples, we have to check 
the intrinsic luminosity distribution of each type of Seyfert 
galaxy using some isotropic radiation. Here we investigate 
the distributions of the mid- and far-infrared luminosity, 
i.e., IRAS 25 /im and 60 /im luminosities. 1 These luminosi- 
ties are thought to scale the nuclear continuum luminos- 
ity which is absorbed and re-radiated by dusty tori, and 
to have little viewing-angle dependence (e.g., Pier, Kro- 
lik 1992; Efstathiou, Rowan-Robinson 1995; Fadda et al. 
1998). Note that the effect of star formation at 25 /zm is 
less than at 60 /zm though active starburst may contribute 
to the 25 /im luminosity, not only the 60 /im luminosity. 
The histograms of the frequency distribution, the median, 
the average, and the 1 a deviation of the IRAS 25 /im and 
60 /im luminosities are given in figures 2 and 3, and table 2. 
The KS test leads to the result that the samples of Sis and 
SI. 5s have higher IRAS 25 /im luminosity than that of the 
S2s, though the samples are statistically indistinguishable 
in the IRAS 60 /im luminosity. Since the intrinsic luminos- 
ity of Seyfert nuclei is more accurately represented by the 
IRAS 25 /im luminosity, this may mean that the samples 
of Sis and SI. 5s are more luminous than that of S2s, sta- 
tistically. We check whether or not this difference in the 
intrinsic luminosity distribution affect the later discussions 
in subsection 3.2. 

3. RESULTS 
3.1. Emission-Line Ratios with Balmer lines 

Historically, the diagnostic diagrams proposed by 
Vcillcux and Osterbrock (1987; hereafter V087) have often 
been used to examine the physical properties of gas clouds 
in NLRs. The V087 diagrams are made from the emission- 
line flux ratios of a forbidden line to a narrow compo- 
nent of a hydrogen Balmer line; e.g., [O ni]A5007/H/3, [N 
n]A6583/Ha, and so on. However, since it is often difficult 
to measure the narrow Balmer components for Sis (and 
SI. 5s) accurately, it is unclear whether or not these flux 
ratios of Sis (and SI. 5s) can be used to study the prop- 
erties of NLRs and can be compared with those of S2s. 
Therefore, prior to comparing the frequency distributions 
of these emission-line flux ratios among various types of 
Seyfert galaxies, we examine whether or not the V087 di- 
agrams can work even for Sis and SI. 5s. 

In figure 4, we show the diagnostic diagrams proposed 
by V087, i.e., the diagrams of the emission- line flux ra- 
tio of [O iii]A5007/H/3 versus that of [N n]A6583/Ha, [S 
n]AA6717,6731/Ha, and [O i]A6300/Ha. The compiled 
data of Seyfert galaxies are plotted in these diagrams with 
the data of extragalactic H n systems for references. Note 

is assumed. The data of the IRAS 25 fim and 60 (im are taken from 
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that these emission-line ratios are little influenced by dust 
extinction, because the wavelength separations between 
the concerned two lines are small. It appears that the data 
of the Sis and the SI. 5s show larger scatters than those 
of the S2s in each diagram, especially in the diagram of 
[O Hl]A5007/H/3 versus [S ii]A6717,6731/Ha (figure 4b). 
Does this suggest that there is a systematic difference in 
the physical properties of NLR gas clouds between the S2s 
and the others, or that the measurements of the narrow 
components of Balmer lines are not well determined for 
the Sis and the SI. 5s? To examine this issue, we inves- 
tigate the relationship of the emission-line flux ratio be- 
tween H7/H/3 and Hcu/H/3 (figure 5). When we assume the 
case B approximation, their theoretically predicted ratios 
are Ha/H/3 = 2.9 and H 7 /H/3 = 0.47 (see, e.g., Ostcr- 
brock 1989). 2 In figure 5, although most of the S2s can 
be well described by the case B prediction taking the ef- 
fects of dust extinction into account (approximately mag 
<; Ay <J 3 mag), most of the Sis cannot be explained by 
the case B prediction with dust extinction. This suggests 
that the fluxes of the narrow components of Balmer lines 
are not well determined for the Sis (and SI. 5s) although 
other possibilities (e.g., contribution of optically-thin gas 
clouds) cannot be ruled out. In any case, it is safe to 
avoid the fluxes of the narrow components of the Balmer 
lines in order to investigate the physical properties of NLR 
gas clouds in Sis and in SI. 5s and in order to study any 
systematic differences of the NLRs between Sis and S2s. 

3.2. Emission-Line Ratios without Balmer Lines 

Some other diagnostic diagrams in which the narrow 
Balmer emission is not used have been proposed to dis- 
cuss the physical properties of gas clouds in Scyfert nuclei 
(e.g., Baldwin et al. 1981; Ohyama 1996; Nagao et al. 
2001a; see also Halpern, Stcincr 1983). Such diagnostic 
diagrams seem to be highly useful in investigating both 
the physical properties of NLRs of Sis and SI. 5s and any 
systematic difference from those of S2s. Therefore, we 
consider statistical properties of the NLRs only by using 
various forbidden emission-line flux ratios (see also Nagao 
et al. 2001a). 

In figures 6a-p, we show the frequency distributions of 
the compiled emission- line ratios, [O i]A6300/[O m]A5007, 
[O n]A3727/[0 m]A5007, [O i]A6300/[O n]A3727, [O 
m]A4363/[0 m]A5007, [S n]A6717/[S n]A6731, [O 
i]A6300/[S n]AA6717,6731, [O n]A3727/[S n]AA6717,6731, 
[S n]AA6717,6731/[0 m]A5007, [O i]A6300/[N n]A6583, 
[O n]A3727/[N n]A6583, [N n]A6583/[0 m]A5007, [S 
n]AA6717,6731/[N n]A6583, [Ne m]A3869/[0 m]A5007, 
[Ne m]A3869/[0 n]A3727, [Ne v]A3426/[0 n]A3727, and 
[Fe vn]A6087/[O m]A5007 for the Sis, the S1.5s and the 
S2s. In these figures, the distributions of the flux ratios 
for the NLSls, the BLSls, the S2 RBL rs, the S2 HBL rs, 
and the S2 _ s are also shown in the right-hand panels. 
The frequency distributions of [O n]A7325/[0 n]A3727, 
[S m]A9069/[S n]AA6717,6731, [N i]A5199/[N n]A6583, 
[N n]A5755/[N n]A6583, [Ar m]A7136/[0 m]A5007, [Fc 
x]A6374/[0 m]A5007, and [Fc xi]A7892/[0 m]A5007 for 
the Sis, the SI. 5s and the S2s are shown in figure 7. Since 

2 These predictions are insensitive to gas temperature and hydrogen density except for a very high-dense gas clouds as existing in BLRs. In 
the range of 5000 K < T < 20000 K and 10 2 cm~ 3 < n H < 10 4 cm~ 3 , these predicted flux ratios vary within 10 % (see Osterbrock 1989). 
Note that, in NLRs of Scyfert galaxies, the Ha emission may be enhanced by collisional excitation (e.g., Ferland, Netzer 1983; Halpern, Stcincr 
1983). 



these line ratios have been measured for a subset of the 
samples, we do not show their histograms for each sub- 
class of Sis and S2s in this figure. In table 3, the median, 
the average, and the la deviation of each emission-line flux 
ratio for the Sis, the Si. 5s and the S2s are given. 

In order to check whether or not the frequency dis- 
tributions of these emission-line flux ratios are statisti- 
cally different among the Sis, the SI. 5s and the S2s, 
we apply the KS test. The resultant KS probabilities 
are given in table 4. These results can be summa- 
rized as follows. (1) As for the emission-line flux ra- 
tios of [O i]A6300/[O m]A5007, [O n]A3727/[0 m]A5007, 
[O i]A6300/[O n]A3727, [S n]A6717/[S n]A6731, [S 
m]A9069/[S n]AA6717,6731, [O i]A6300/[S n]AA6717,6731, 
[O n]A3717/[S n]AA6717,6731, [N i]A5199/[N n]A6583, 
[O i]A6300/[N n]A6583, [O n]A3727/[N n]A6583, [N 
n]A6583/[0 m]A5007, [S n]AA6717,6731/[N n]A6583 and 
[Ar m]A7136/[0 m]A5007, there are little or no sys- 
tematic differences among the Sis, the SI. 5s and the 
S2s. (2) The emission-line flux ratios of [O n]A7325/[0 
n]A3727, [N n]A5755/[N n]A6583 and [Fe xi]A7892/[0 
iii]A5007 appear to be slightly higher in the Sis than 
those in the S2s, although these differences are statisti- 
cally insignificant. (3) It is statistically significant that 
the emission- line flux ratios of [O iii]A4363/[0 iii]A5007, 
[S ii]AA6717,6731/[0 iii]A5007, [Ne m]A3869/[0 m]A5007, 
[Ne m]A3869/[0 ii]A3727, [Ne v]A3426/[0 n]A3727, [Fe 
vn]A6087/[O m]A5007 and [Fe x]A6374/[0 m]A5007 are 
different between the S2s and the other types of Seyfcrt 
galaxies (i.e., the Sis and the SI. 5s). Since these 23 
emission-line flux ratios are not correlated with the red- 
shift and IRAS 25 /im luminosity as shown in figure 8, the 
above results seem to be almost free from the redshift and 
luminosity biases mentioned in subsection 2.2. 

As mentioned in subsection 2.1, no reddening correc- 
tion has been made for all of the compiled emission-line 
flux ratios. Since it is known that the dust extinction is 
larger on average in S2s than that in Sis (e.g., Dahari, De 
Robertis 1988), the above results may be caused by the dif- 
ference in the amount of extinction between the Sis and 
the S2s. In order to check whether or not this is the case, 
we examine the effect of the extinction correction for each 
emission-line flux ratio adopting the Cardelli's extinction 
curve (Cardelli et al. 1989). In table 5, we summarize the 
correction factors, by which the emission-line flux ratios 
should be multiplied to be converted into the extinction- 
corrected values for the case of Ay = 1 mag. Note that the 
mean difference of the amount of dust extinction between 
Sis and S2s is ~1 mag (Dahari, Dc Robertis 1988; see also 
Dc Zotti, Gaskell 1985). Since the effect of the extinction 
correction is too small for the cases of [O m]A4363/[0 
m]A5007, [Ne m]A3869/[0 m]A5007 and [Ne v]A3426/[0 
n]A3727, the differences in these emission-line flux ratios 
between the Sis and the S2s cannot be attributed only to 
the effect of the dust extinction. Furthermore, the differ- 
ences in the flux ratios of [Ne m]A3869/[0 n]A3727, [Fe 
vn]A6087/[O ni]A5007 and [Fe x]A6374/[0 m]A5007 can- 
not be also interpreted by the effect of the dust extinction. 
Therefore, we conclude that the AGN-type dependence of 
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these emission-line flux ratios is due not to the difference 
in the amounts of extinction, but to some other factors, as 
discussed later. 

Here, we mention that these results are consistent with 
the previous studies. The excess of the flux ratio of [O 
lll]A4363/[0 iii]A5007 in Sis (and in SI. 5s) has been re- 
ported by Osterbrock et al. (1976), Heckman and Balick 
(1979), Shuder and Osterbrock (1981), Cohen (1983), and 
Nagao et al. (2001b). The excess of the intensities of high- 
ionization iron emission lines of Sis has been mentioned 
by Shuder and Osterbrock (1981), Cohen (1983), Mu- 
rayama and Taniguchi (1998a), and Nagao et al. (2000). 
Schmitt (1998) noted that Sis exhibit higher ratios of [Nc 
m]A3869/[0 n]A3727 and [Ne v]A3426/[0 n]A3727 than 
S2s. 

4. DISCUSSION 

4.1. The Origin of the Seyfert-Type Dependence of the 
Emission-Line Flux Ratios 

The results presented in the last section can be sum- 
marized as follows: (1) Most of the emission-line flux ra- 
tios which show the Seyfert-type dependence contain a 
high critical-density (p> 10 7 cm -3 ) and/or high ionization- 
potential (> 100 eV) emission lines, except for the ratio 
of [S Il]AA6717,6731/[0 m]A5007. (2) On the other hand, 
the flux ratios which do not exhibit the Seyfert-type de- 
pendence do not contain high-ionization emission lines. In 
table 6, we summarize the ionization potential and the 
critical density of each emission line used in our analysis. 
Since the [Fe xi]A7892 emission has a very high ionization 
potential, its relative intensity could be different between 
Sis and S2s. However, it is unclear in our sample whether 
or not the [Fe xi]A7892 emission is stronger in the Sis than 
in the S2s. One of the reasons for this may be the small 
number of [Fe Xl] -detected objects. Further observations 
will be necessary to confirm the difference in the frequency 
distribution of the [Fe xi]A7892 intensity between Sis and 
S2s. 

We now consider the origin of the Seyfert-type depen- 
dence of those emission- line flux ratios. There are two 
possible alternatives. One is that the high-ionization emis- 
sion lines arise mainly from dense gas clouds which are 
located very close to nuclei (Torus HINER 3 ; see Mu- 
rayama, Taniguchi 1998a, b). Since such a region can 
be hidden by dusty tori, the visibility of the dense gas 
clouds may depend on a viewing angle toward the tori. 
This component may correspond to highly-ionized, dense 
(~ 10 7 ~ 8 cm~ 3 ) gas beside the inner wall of dusty tori 
(see Pier, Voit 1995). This idea was proposed by Mu- 
rayama and Taniguchi (1998a) to explain the stronger 
[Fe Vii]A6087 emission in Sis compared to S2s. Taking 
this Torus-HINER component into account, Murayama 
and Taniguchi (1998b) showed that the difference in the 
[Fe vn]A6087 intensity between Sis and S2s can be suc- 
cessfully interpreted by a dual-component photoionization 
model proposed by them (see also Nagao et al. 2001b). 
The other idea is that the difference in the emission-line 
flux ratios reflects the intrinsic difference of NLR proper- 
ties, e.g., physical size, density, temperature, and so on. 
Osterbrock (1978) mentioned that the systematic differ- 
ence in the flux ratio of [O m]A4363/[0 m]A5007 can 

3 The term "HINER" means high-ionization nuclear emission-line 



be understood assuming njj ~ 10 6 ~ 7 cm~ 3 for Sis and 
njj < 10 5 cm~ 3 for S2s. On the contrary, Heckman and 
Balick (1979) and Cohen (1983) claimed that the origin 
of the difference in [O m]A4363/[0 m]A5007 between Sis 
and S2s is attributed not to the density difference, but 
to the temperature difference; i.e., T c > 2 x 10 4 K for 
Sis while T c ~ 10 4 K for S2s. Here, it must be noted 
that these situations may occur when high-density or high- 
temperature regions are hidden by any obscuring matter 
in S2s. Namely, these two scenarios do not necessarily 
mean the intrinsic difference in the NLR properties. On 
the other hand, Schmitt (1998) reported that the system- 
atic differences in the ratios of [Ne m]A3869/[0 n]A3727 
and [Ne v]A3426/[0 n]A3727 between Sis and S2s can 
be explained by taking the physically (i.e., not projected) 
smaller NLR size of Sis compared to that of S2s, as sug- 
gested by Schmitt and Kinney (1996), into account. Since 
smaller NLRs contain fewer ionization-bounded clouds, 
which radiate low-ionization emission lines selectively, this 
results in weaker [O n]A3727 emission compared to [Ne 
lll]A3869 and [Ne v]A3426 in Sis (hereafter "smaller NLR 
model" ) . 

Which scheme is more realistic, the obscuration of the 
Torus HINER in S2s or the smaller NLR model? In either 
case, both high-ionization gas clouds and low-ionization 
ones are necessary to explain the observations, and a cer- 
tain systematic difference in the relative contribution be- 
tween the two components can be regarded as the origin 
of the Seyfert-type dependence of the emission-line flux 
ratios. However, the two ideas make different predictions. 
The Torus-HINER model predicts that the intrinsic lumi- 
nosity of high-ionization lines are brighter in Sis than in 
S2s because the Torus HINER is assumed to be hidden in 
the S2s (Murayama, Taniguchi 1998a). On the other hand, 
the smaller NLR model predicts that the intrinsic lumi- 
nosities of low-ionization lines are lower in the Sis than in 
the S2s because the low-ionization lines arise mostly from 
ionization-bounded gas clouds at outer NLRs. Therefore, 
the Torus-HINER model predicts similar [O m]A5007 lu- 
minosity among the Scyfert types and higher luminosities 
of [O m]A4363, [Ne m]A3869, and [Fe vn]A6087 in Sis than 
in S2s, while the "smaller NLR model" predicts lower [O 
ni]A5007 luminosity in Sis than that in S2s, and similar lu- 
minosities of [O m]A4363, [Ne m]A3869, and [Fe vn]A6087 
among the Seyfert types. 

We now discuss observational tests for the two models. 
Since the frequency distribution of the intrinsic luminosity 
is different among the Sis, the SI. 5s and S2s in our sample, 
as mentioned in subsection 2.2, these emission-line lumi- 
nosities should be normalized by IRAS 25 fira luminosity 
to test the above predictions. In figure 9, we show the fre- 
quency distributions of the [O m]A5007, [O m]A4363, [Ne 
m]A3869, and [Fe vn]A6087 luminosities of the Sis, the 
SI. 5s, and the S2s, which are normalized by the IRAS 25 
/im luminosity. The median, the average, and the la devi- 
ation of each relative luminosity for the Sis, the SI. 5s, and 
the S2s are given in table 7. In order to examine whether 
or not the distribution of the relative strength of the emis- 
sion lines depends on the Seyfert type, we apply the KS 
test. The resultant KS probabilities are given in table 8. 
The KS test leads to the following results: 

region (Binettc 1985; Murayama et al. 1998) 



Nagao, Murayama, & Taniguchi 



• The [O m]A5007 luminosity normalized by IRAS 
25 /im luminosity is statistically indistinguishable 
among the Seyfert types. This is consistent with 
the Torus-HINER model, but in conflict with the 
prediction of the "smaller NLR model" . 

• The [O Hl]A4363 and the [Fe Vii]A6087 luminosity 
normalized by IRAS 25 /im luminosity appear to be 
higher in the Sis and in the SI. 5s than in the S2s, al- 
though the statistical significance is low. This agrees 
with the prediction of the Torus-HINER model, but 
cannot be understood in terms of the "smaller NLR 
model" . 

• The [Ne m]A3869 luminosity normalized by IRAS 
25 /xm luminosity is statistically indistinguishable 
among the Seyfert types. This seems to disagree 
with the prediction of the Torus-HINER model. 

The first two results seem to support the Torus-HINER 
model. As for the last result, we have to take into account 
that the critical density of the [Ne iii]A3869 transition is 
smaller than those of the [O m]A4363 and [Fe vn]A6087 
transitions. This may be the reason why the frequency 
distribution of the relative luminosity of the [Ne in] A3869 
emission is indistinguishable among the Seyfert types. We 
therefore conclude that the Torus-HINER model appears 
to be much more consistent with the observation than the 
smaller NLR model. We thus adopt the Torus-HINER 
model in the following discussion. 

One of our interesting results is that the emission-line 
flux ratio of [S n]AA6717,6731/[0 m]A5007 exhibits a sys- 
tematic difference in its frequency distribution between the 
SI. 5s and the S2s. If this difference is also attributed to dif- 
ferent contributions of the Torus HINER between the two 
types, it is suggested that the visibility of the [O iii]A5007 
emitting region depends on the viewing angle toward dusty 
tori. Indeed, the isotropy of the [O iii]A5007 emission has 
sometimes been called into question, especially in radio- 
loud AGNs. Jackson and Browne (1990) reported that the 
[O iii]A5007 luminosity of the narrow- line radio galaxies is 
lower by 5-10 times than that of the broad-line quasars, 
matched in redshift and extended radio luminosity. On 
the other hand, Hes et al. (1993) found that these two 
kinds of radio-loud AGNs show no difference in the [O 
n]A3727 luminosity. These results suggest that the [O 
n]A3727 emission is isotropic, but the [O iii]A5007 emis- 
sion is not (see also Baker, Hunstead 1995; Baker 1997; 
cf., Simpson 1998). Note that this conclusion is based 
on an assumption that broad-line quasars and (powerful) 
narrow-line radio galaxies are intrinsically similar, but dif- 
ferent only in the viewing angle (see, e.g., Barthcl 1989). 
Polarized [O iii]A5007 emission has been detected in some 
radio galaxies (di Scrego Alighicri et al. 1997), supporting 
the scenario that a part of the [O iii]A5007 flux is hidden 
by the tori. Because such polarized [O iii]A5007 emission 
has also been detected in a Seyfert 2 galaxy NGC 4258 
(Wilkes et al. 1995; Barth et al. 1999), it is interesting to 
investigate whether or not the [O iii]A5007 emission has an 
anisotropic property also in the sample of Seyfert galaxies, 
although Mulchaey et al. (1994) reported a negative result 
for this possibility. 

If the [O iii]A5007 emission of our sample is affected by 
an orientation-dependent dust obscuration, the frequency 



distribution of L([0 iii]A5007) would show a larger scatter 
than that of L([0 Il]A3727) as a result of the obscura- 
tion of [O iii]A5007 emission at a large inclination angle 
(sec Kuraszkiewicz et al. 2000). Here, the [O n]A3727 
emission is thought to have no viewing-angle dependence. 
In figure 9, we show a diagram of L([0 iii]A5007) versus 
L([0 n]A3727), in which no obscuration effect is found. 
Recently, Kuraszkiewicz et al. (2000) reported that radio- 
quiet quasars do not show evidence for an anisotropic 
property of the [O lll]A5007 emission. Therefore, tak- 
ing the results of both this study and Kuraszkiewicz 
et al. (2000) into account, radio-quiet AGNs including 
quasars and Seyfert galaxies may not have the anisotropic 
property of the [O m]A5007 emission, contrary to radio- 
loud AGNs (see also Mulchaey et al. 1994). Then, 
why is the frequency distribution of the flux ratio of [S 
n]AA6717,6731/[0 m]A5007 different between the SI. 5s 
and the S2s? A possible reason may be contamination 
of lower-ionization emission-line fluxes arising from cir- 
cumnuclear star-forming regions associated in S2s, be- 
cause it is known that circumnuclear star-forming activity 
tends to associate with S2s more frequently than with Sis 
(e.g., Thuan 1984; Heckman et al. 1995, 1997). Indeed, 
there is a weak negative correlation between the flux ra- 
tio of [S n]AA6717,6731/[0 m]A5007 and that of vF„(25 
/iTa)/h'F l/ (60 /im), as shown in figure 11. Since the ra- 
tio of isF l/ (25/im)/h'F l ,(60iim) becomes smaller when the 
star-formation activity contributes much more to the in- 
frared continuum radiation, this correlation suggests that 
the star-formation activity may enhance the flux ratio of 
[S n]AA6717,6731/[0 m]A5007. 

4.2. Physical Properties of the Torus HINERs 

Although the Torus HINER is regarded as a "highly- 
ionized and dense" component, there are some remaining 
issues to be investigated: How dense is the Torus HINER? 
And, how high is the ionization parameter of ionized gas 
in the Torus HINER? To examine these issues, we carry 
out single cloud photoionization model calculations using 
the spectral synthesis code Cloudy version 90.04 (Ferland 
1996), which solves the equations of statistical and ther- 
mal equilibrium and produces a self-consistent model of 
the run of temperature as a function of depth into a neb- 
ula. Here, we assume that a uniform density, dust-free gas 
cloud with plane-parallel geometry is ionized by a power- 
law continuum source. The parameters for the calcula- 
tions arc: (1) the spectral energy distribution (SED) of 
the input radiation, (2) the chemical composition, (3) the 
hydrogen density of the cloud («h), (4) the ionization pa- 
rameter (U), and (5) the hydrogen column density. We 
adopt the following input continuum spectrum: (i) a = 
2.5 for A > 10/im, (ii) a = -2.0, -1.5, and -1.0, between 10 
/im and 50 keV, and (iii) a = -2.0 for hv > 50 keV, where 
a is a spectral index (/„ oc v a ). Koski (1978) reported 
that the optical continua of S2s can be approximated by 
a stellar contribution diluted by a featureless continuum, 
with the latter component described by a power law with 
a = —1.5 ±0.5 (see also Storchi-Bergmann, Pastoriza 
1989, 1990; Kinney et al. 1991). We set the gas-phase 
elemental abundances to be the solar ones taken from 
Grevesse and Anders (1989) with extensions by Grevesse 
and Noels (1993). To see the effects of mctallicity, we also 
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calculated for the case of Z — 0.5 and 2.0 in addition to 
Z = 1.0. We performed several model runs covering 10 2 
cm" 3 < n H < 10 8 cm" 3 and 10" 35 < U < 10 -1 ' 5 for three 
kinds of SED. Since an unusually strong [O i]A6300 emis- 
sion with respect to the [O III]A5007 emission is predicted 
by the model assuming high-density, ionization-bounded 
gas clouds, we assumed truncated (i.e., matter-bounded) 
clouds (see, e.g., Murayama, Taniguchi 1998b). To make 
the gas clouds matter-bounded ones, the calculations were 
stopped at a hydrogen column density when a Lyman limit 
optical depth (rg^) reached up to 0.1. In this condition, 
more than 96 % of hydrogen is ionized at the most outer 
region of a cloud. 

In figure 12, the calculated emission- line flux ratios of [O 
m]A4363/[0 m]A5007, [Ne m]A3869/[0 m]A5007 and [Fe 
vil]A6087/[O iii]A5007 are plotted as a function of n H . In 
order to explain the observed flux ratio of [O m]A4363/[0 
lll]A5007 (= 0.122 ± 0.116 for the Sis), rather high density 
(riH ~ 10 6 cm" 3 ) is necessary, regardless of U, a, and the 
metallicity. To examine the ionization parameter which 
can consistently explain the observed emission-line flux 
ratios of [Ne m]A3869/[0 m]A5007 and [Fe vn]A6087/[O 
m]A5007 (0.246 ± 0.145 and 0.122 ± 0.116 for the Sis, 
respectively) adopting rtn ~ 10 6 cm" 3 , we show the cal- 
culated flux ratios of [Ne m]A3869/[0 m]A5007 and [Fe 
Vll]A6087/[O iii]A5007 as a function of U for the case of 
riH ~ 10 6 cm" 3 in figure 13. This figure suggests that the 
observed emission-line flux ratios can be described by the 
models with U = 10" 2 . This value is almost independent 
of the input SEDs and metallicities. 

However, these estimates (nn ~ 10 6 cm" 3 and U = 
10 -2 ) are based on the assumption that all of the ob- 
served emission-line fluxes arise from the Torus HINER. 
In fact, it is evident that a part of the emission-line flux 
arises from the outer low-density regions. Since the gas 
clouds in such regions are expected to emit relatively low- 
ionization emission-line spectra compared to those in the 
Torus HINER, the Torus-HINER emission must be diluted 
by such low-ionization emission-line spectra. Therefore, 
we should remind that the derived properties of Torus 
HINERs are lower limits, i.e., nn > 10 6 cm" 3 and 
U > 10" 2 . 

These results are not modified significantly if we take a 
smaller value of T912. In figure 14, we show the calculated 
emission- line flux ratios of [O m]A4363/[0 m]A5007, [Ne 
m]A3869/[0 m]A5007, and [Fe vn]A6087/[O m]A5007 in 
the case of T912 = 0.01 as a function of tih (only the mod- 
els with solar abundances are shown). The results of the 
calculations are nearly the same as those for the case of 
T912 =0.1. On the contrary, relatively smaller flux ratios 
of [Fe vn]A6087/[O m]A5007 are predicted if we take a 
larger value of T912 . For reference we show the calculation 
results for the case of r 9 i 2 = 1.0 in figure 15. Accordingly, 
larger gas density and/or ionization parameters than those 
estimated for the case of T912 =0.1 are necessary to ex- 
plain the observations in the cases of Tgi2 > 0.1. In any 
case, we conclude that the lower limits of the gas den- 
sity and the ionization parameter of the Torus-HINER are 
10 6 cm" 3 and U ~ 10" 2 . 



«H 



5. SUMMARY 



Based on a compilation of the optical emission-line spec- 
tra of Seyfert galaxies, we have investigated the Seyfert- 



type dependences of various emission-line flux ratios. Our 
analysis was made using only forbidden emission lines. 
This method has enabled us to compare the physical prop- 
erties of the NLR among the various types of Seyfert nuclei 
(e.g., SI, SI. 5, and S2). 

In consequence of the statistical comparisons of vari- 
ous forbidden emission-line flux ratios among the Seyfert 
types, we have obtained the following results: 

• As for the emission-line flux ratios of [O i]A6300/[O 
m]A5007, [O n]A3727/[0 m]A5007, [O i]A6300/[O 
n]A3727, [S n]A6717/[S n]A6731, [S m]A9069/[S 
n]AA6717,6731, [O i]A6300/[S n]AA6717,6731, [O 
n]A3717/[S n]AA6717,6731, [N i]A5199/[N n]A6583, 
[O i]A6300/[N n]A6583, [O n]A3727/[N n]A6583, [N 
n]A6583/[0 m]A5007, [S n]AA6717,6731/[N n]A6583 
and [Ar m]A7136/[0 m]A5007, there are little or no 
systematic differences among the Sis, the SI. 5s and 
the S2s. 

• On the other hand, it is statistically significant 
that the emission- line flux ratios of [O iii]A4363/[0 
m]A5007, [S n]AA6717,6731/[0 m]A5007 [Ne 
m]A3869/[0 m]A5007, [Ne m]A3869/[0 n]A3727, 
[Ne v]A3426/[0 n]A3727, [Fe vn]A6087/[O m]A5007 
and [Fe x]A6374/[0 m]A5007 are systematically 
higher in the Sis and the SI. 5s than in the S2s. 

These results can be interpreted as that the flux ratios of 
a high-ionization and a high-critical-density emission line 
to a low-ionization emission line are systematically higher 
in the Sis than in S2s. 

There are two ideas which can possibly explain the 
Seyfert-type dependences of the relative strength of high- 
ionization emission lines; i.e., the viewing-angle depen- 
dence of visibility of Torus HINER and the intrinsic dif- 
ference of the NLR size between Sis and S2s. To dis- 
criminate these two possibilities, we have examined the 
Seyfert-type dependences of some emission-line luminosi- 
ties. The results are consistent with the prediction of the 
Torus-HINER model, but disagree with the "smaller NLR 
model". We have concluded that the difference of visi- 
bility of Torus HINER between Sis and S2s causes the 
Seyfert-type dependences of the relative strength of the 
high-ionization emission lines. 

In order to investigate the properties of the Torus 
HINER, we have compared the relative strengths of high- 
ionization emission lines with the results of the photoion- 
ization model calculations. The estimated lower limit of 
the density and the ionization parameter are «h ~ 10 6 
cm" 3 and U ~ 10" 2 . These constraints are almost in- 
dependent of input SEDs and metallicities of the ionized 
gas. 
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Table 1. Medians, averages, and 1 a deviations of redshift, vL v (2h /im), and i , L l/ (60 /im) for each type of Seyfert 

galaxies. 







Redshift 


log 


vL v {2h /an) 


log 


vL„(60 /im) 


Type 


Median 


Average and 1 a 


Median 


Average and 1 a 


Median 


Average and 1 a 


Sltotal 


0.0502 


0.1023 ± 0.1306 


42.296 


42.315 ± 0.555 


42.179 


42.203 ± 0.629 


NLS1 .. 


0.0532 


0.0687 ± 0.0633 


42.207 


42.166 ± 0.672 


42.065 


42.235 ± 0.632 


BLS1 .. 


0.0500 


0.1138 ± 0.1452 


42.297 


42.383 ± 0.490 


42.189 


42.190 ± 0.636 


SI. 5 


0.0384 


0.0985 ± 0.1514 


42.553 


42.339 ± 0.639 


42.240 


42.231 ± 0.606 


S2 to tal • • • • 


0.0202 


0.0373 ± 0.0660 


41.933 


41.795 ± 0.765 


42.000 


41.979 ± 0.653 


S2rblr 


0.0133 


0.0437 ± 0.0967 


41.545 


41.432 ± 0.820 


41.735 


41.728 ± 0.617 


S2hblr 


0.0136 


0.0214 ± 0.0185 


42.430 


42.447 ± 0.322 


42.176 


42.213 ± 0.523 


S2" ... 


0.0240 


0.0364 ± 0.0509 


41.933 


41.883 ± 0.693 


42.093 


42.061 ± 0.661 



Table 2. Resultant KS probabilities* concerning redshift, vL u (25 /im), and vL v (Q0 /im) 



Type Sl tot ai NLS1 BLS1 


SI. 5 S2total 


S2rblr 


S2hblr 


S2~ 




Redshift 


Sltotal 






1.232XKT 1 1.291xl(T 12 


7.673xl(T lu 


5.925xl(T 4 


2.343 xl(T' 


NLS1 




6.382XKT 1 


6.153XKT 1 4.563xl(T 5 


1.237xl(T 5 


3.103xl(T 3 


1.688xl(T 3 


BLS1 






1.466XKT 1 4.865XKT 11 


2.515xl(T 9 


6.994xl(T 4 


1.068xl(T 6 


S1.5 






1.458xl(T 7 


5.421xl(T 7 


8.856xl(T 3 


1.248xl(T 4 


S2 to taI 














S2rblr 










5.264XKT 1 


1.656xl(T 2 


S2hblr 












8.871xl(T 2 


S2" 


















vLv(25 /im) 


Sltotal • 






6.952xl(T 2 5.904xl(T 4 


1.816xl(T 4 


5.321XKT 1 


2.534xl(T 3 


NLS1 . 




3.223X1CT 1 


2.371X1CT 1 1.724XKT 1 


2.475xl(T 2 


3.242XKT 1 


2.987XKT 1 


BLS1 .. 






1.238XKT 1 1.679xl(T 4 


4.898 xl(T 


6.685XKT 1 


6.470 xl(T 4 


S1.5 ... 






4.121xl(T 4 


2.000xl(T 4 


7.709XKT 1 


1.483 xl(T 3 


S2 to tal • 














S2rblr 










1.184xl(T 3 


1.642XKT 1 


S2hblr 












1.037xl(T 2 


S2- ... 


















vL v (6Q /im) 


Sltotal • 






7.619XKT 1 5.383xl(T 2 


2.389xl(T 3 


9.811XKT 1 


3.852XKT 1 


NLS1 . 




8.215XKT 1 


9.358XKT 1 6.148XKT 1 


1.230XKT 1 


9.585XKT 1 


9.072XKT 1 


BLS1 .. 






5.805XKT 1 2.455xl(T 2 


1.834xl(T 3 


9.958XKT 1 


1.911xl0 _1 


S1.5 ... 






1.044XKT 1 


9.109xl(T 3 


8.899XKT 1 


4.809XKT 1 


S2 to tal • 














S2rblr 










1.252XKT 1 


8.458xl(T 2 


S2hblr 












7.279XKT 1 


S2- ... 















These values means the probabilities of two sample drawn from the same parent population. 
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Table 3. Median, average, and 1 a deviation of each line ratio. 



Line Ratio 


Seyfert 1 Galaxies 


Seyfert 


1.5 Galaxies 


Seyfert 2 Galaxies 




Median 


Average and 1 a 


Median 


Average and 1 a 


Median 


Average and 1 a 


[0 i]A6300/[O m]A5007 


0.074 


0.129 ± 0.175 


0.064 


0.110 ± 0.165 


0.091 


0.135 ± 0.149 


[0 n]A3727/[0 m]A5007 


0.201 


0.280 ± 0.200 


0.165 


0.223 ± 0.206 


0.233 


0.341 ± 0.396 


[0 i] A6300/[O m]A3727 


0.322 


0.447 ± 0.378 


0.397 


0.522 ± 0.402 


0.331 


0.514 ± 1.048 


[0 m]A4363/[0 m]A5007 .... 


0.073 


0.101 ± 0.101 


0.063 


0.104 ± 0.200 


0.025 


0.038 ± 0.043 


[0 n]A7325/[0 n]A3727 


0.340 


0.496 ± 0.639 


0.195 


0.231 ± 0.099 


0.239 


0.244 ± 0.112 


[S n]A6717/[S n]A6731 


1.025 


1.083 ± 0.307 


1.000 


1.001 ± 0.167 


1.072 


1.086 ± 0.187 


[S m]A9069/[S n]AA6717,6731 


0.692 


0.617 ± 0.315 


0.481 


0.538 ± 0.209 


0.567 


0.605 ± 0.208 


[0 i]A6300/[S n]AA6717,6731 


0.302 


0.446 ± 0.392 


0.359 


0.483 ± 0.352 


0.267 


0.333 ± 0.274 


[0 n]A3727/[S n]AA6717,6731 


0.903 


1.211 ± 1.103 


0.891 


0.990 ± 0.545 


0.853 


0.947 ± 0.485 


[S n]AA6717,6731/[0 m]A5007 


0.209 


0.330 ± 0.285 


0.180 


0.242 ± 0.238 


0.332 


0.555 ± 0.704 


[N i]A5199/[N n]A6583 


0.037 


0.045 ± 0.025 


0.035 


0.051 ± 0.047 


0.030 


0.033 ± 0.017 


[N n]A5755/[N n]A6583 


0.028 


0.048 ± 0.055 


0.022 


0.027 ± 0.017 


0.012 


0.013 ± 0.009 


[0 i]A6300/[N n]A6583 


0.166 


0.264 ± 0.303 


0.271 


0.395 ± 0.725 


0.136 


0.180 ± 0.126 


[0 n]A3727/[N n]A6583 


0.432 


0.635 ± 0.622 


0.454 


0.550 ± 0.372 


0.424 


0.593 ± 0.459 


[N n]A6583[0 m]A5007 


0.459 


0.820 ± 0.967 


0.311 


0.486 ± 0.554 


0.570 


1.080 ± 1.523 


[S n]AA6717,6731/[N n]A6583 


0.462 


0.724 ± 0.665 


0.550 


0.708 ± 0.472 


0.541 


0.613 ± 0.450 


[Ar iii]A7136/[0 iii]A5007 . . . 


0.029 


0.036 ± 0.021 


0.021 


0.029 ± 0.021 


0.041 


0.054 ± 0.032 


[Ne m]A3869/[0 m]A5007 . . . 


0.210 


0.246 ± 0.145 


0.134 


0.152 ± 0.088 


0.080 


0.110 ± 0.118 


[Ne m]A3869/[0 n]A3727 .... 


0.965 


1.178 ± 0.951 


0.674 


0.873 ± 0.556 


0.358 


0.449 ± 0.327 


[Ne v]A3426/[0 n]A3727 


1.280 


1.467 ± 1.031 


1.006 


1.519 ± 1.373 


0.370 


0.472 ± 0.433 


[Fe vn]A6087/[O m]A5007 . . . 


0.089 


0.122 ± 0.116 


0.047 


0.082 ± 0.092 


0.017 


0.022 ± 0.020 


[Fe x]A6374/[0 m]A5007 .... 


0.070 


0.077 ± 0.061 


0.035 


0.045 ± 0.034 


0.009 


0.021 ± 0.034 


[Fe xi]A7892/[0 m]A5007 .... 


0.067 


0.073 ± 0.053 


0.026 


0.049 ± 0.057 


0.002 


0.021 ± 0.036 
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Table 4. Resultant KS probabilities* concerning the forbidden emission-line flux ratios. 



SI. 5s vs. S2 



total 



Line Ratio Sl to tai vs. SI. 5s Sl to tai vs. S2 to tai 

[O i]A6300/[O m]A5007 2.066X10- 1 3.140X10" 1 

[O n]A3727/[0 m]A5007 6.928xl(T 2 8. 763x10-* 

[O i]A6300/[O n]A3727 1.467X10- 1 6.653X1CT 1 

[Oiii]A4363/[Oiii]A5007 .... 3.120X10" 1 8.701xl0~ 13 

[O n]A7325/[0 n]A3727 2.006X1CT 1 8.458xl0~ 2 

[S n]A6717/[S n]A6731 2.660X10" 1 5.359X10" 1 

[S m]A9069/[S n]AA6717,6731 3.009X10" 1 4.673X10- 1 

[O i]A6300/[S n]AA6717,6731 3.423X1CT 1 6.035xl0~ 2 

[O n]A3727/[S n]AA6717,6731 6.566X10- 1 3.909X10- 1 

[S n]AA6717,6731/[0 m]A5007 9.093xl(T 2 8.885xl0~ 3 

[N i]A5199/[N n]A6583 8. 958x10^ 4.265X10- 1 

[N n]A5755/[N n]A6583 8. 096x10-* 7.785xl0~ 2 

[O i]A6300/[N n]A6583 1. 588x10^ 2.394X1CT 1 

[O n]A3727/[N n]A6583 8.310X1CT 1 8.702x10-* 

[N n]A6583[0 m]A5007 7.421xl(T 2 5.395x10-* 

[S ii]AA6717,6731/[Nii]A6583 3.472xl0- 2 1.658X1CT 1 

[Ar iii]A7136/[Oiii]A5007 ... 2.826x10-* 2.472x10-* 

[Neiii]A3869/[Oni]A5007 ... 7.646xl(T 5 1.685xl0-* 5 

[Ne m]A3869/[0 n]A3727 .... 9.533xl(T 2 3.494xl0" 13 

[Ne v]A3426/[0 n]A3727 8.275x10-* 2.068xl0- 4 

[Fe vn]A6087/[O m]A5007 ... 1.051x10-* 2.562x10-*° 

[Fcx]A6374/[Om]A5007 .... 2.938xl(T 2 9.535xl0- fi 

[Fe xi]A7892/[Q m]A5007 .... 1.963x10-* 1.006xlCT 2 

These values means the probabilities of two sample drawn from the same parent population. 



2.551xlCT 2 
1.061 xl(T 2 
5.610X10" 2 
2.193xl0- 7 
9.894x10-* 
8.009xl0" 2 
3.168x10-* 
2.670xl0- 3 
8.521x10-* 
5.849xl0~ 5 
4.116x10-* 
1.866xl0- 2 
1.839X10" 3 
9.879x10-* 
2.386xl0~ 3 
4.095x10-* 
6.642xl0- 3 
4.231X10" 7 
1.070xl0" 6 
1.343xl0~ 5 
1.172xl0- 5 
7.356xl0- 4 
1.518xlQ-* 
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Table 5. Correction factors for the extinction (Ay = 1.0 mag) 



Line Ratio Correction Factor 

[O i]A6300/[O m]A5007 0.786 

[O n]A3727/[0 m]A5007 1.471 

[O i]A6300/[O n]A3727 1.872 

[O m]A4363/[0 m]A5007 .... 1.222 

[O n]A7325/[0 n]A3727 0.461 

[S n]A6717/[S n]A6731 1.002 

[S m]A9069/[S n]AA6717,6731 0.745 

[O i]A6300/[S n]AA6717,6731 1.062 

[O n]A3727/[S n]AA6717,6731 1.988 

[S n]AA6717,6731/[0 m]A5007 0.740 

[N i]A5199/[N n]A6583 1.263 

[N n]A5755/[N n]A6583 1.132 

[O i]A6300/[N n]A6583 1.041 

[O n]A3727/[N n]A6583 1.949 

[N ii]A6583[0 m]A5007 0.755 

[S n]AA6717,6731/[N n]A6583 0.980 

[Ar m]A7136/[0 m]A5007 . . . 0.697 

[Ne m]A3869/[0 m]A5007 . . . 1.423 

[Ne m]A3869/[0 n]A3727 .... 0.967 

[Ne v]A3426/[0 n]A3727 1.069 

[Fe vn]A6087/[O m]A5007 . . . 0.811 

[Fc x]A6374/[0 m]A5007 .... 0.778 

[Fe xi]A7892/[Q ni]A5007 . . . . 0.627 
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Table 6. Ionization potentials and the critical densities for forbidden emission lines. 



Line 



Ionization Potential 
Lower (eV) Upper (eV) 



Critical Density 
(cm -3 ) 



[Ne v]A3426 
[O n]A3726 
[O n]A3729 
[Ne m]A3869 
[O m]A4363 
[O m]A5007 
[N i]A5199 
[N n]A5755 
[Fe vn]A6087 
[O i]A6300 
[Fe x]A6374 
[N n]A6583 
[S n]A6717 
[S n]A6731 
[Ar m]A7136 
[O n]A7318.6 
[O n]A7319.9 
[O n]A7329.9 
[O n]A7330.7 
[Fe xi]A7892 
[S iii]A9069 



97.1 
13.6 
13.6 
41.0 
35.1 
35.1 
0.0 
14.5 
99.1 
0.0 
233.6 
14.5 
10.4 
10.4 
27.6 
13.6 
13.6 
13.6 
13.6 
262.1 
23.3 



126.2 
35.1 
35.1 
63.5 
54.9 
54.9 
14.5 
29.6 
125.0 
13.6 
262.1 
29.6 
23.3 
23.3 
40.7 
35.1 
35.1 
35.1 
35.1 
290.1 
34.8 



1.6 xlO 7 

1.4 xlO 4 
3.3 xlO 3 

9.7 xlO 6 

3.3 xlO 7 
7.0 xlO 5 
2.0 xlO 3 
3.2 xlO 7 

3.6 xlO 7 

1.8 xlO 6 

4.8 xlO 9 

8.7 xlO 4 

1.5 xlO 3 

3.9 xlO 3 

4.8 xlO 6 
4.8 xlO 6 

7.4 xlO 6 

8.6 xlO 6 
7.0 xlO 6 

1.2 xlO 6 
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Table 7. Median, average, and 1 a deviation of the line luminosities normalized by h>L u (2§ /zm). 



Seyfert 1 Galaxies Seyfert 1.5 Galaxies Seyfert 2 Galaxies 

Median Average and 1 a Median Average and 1 a Median Average and 1 a 



L([0 in]A5007)/^(25/zm) 
L([0 m]XA363)/vL v (25fj,m) 
L([Ne m]A3869)/i/L J/ (25/xm) 
L([Fe Vil]X6087)/vL v (25(jm) 



0.245 


0.251 ± 0.165 


0.216 


0.536 ± 1.117 


0.179 


0.334 ± 0.403 


0.018 


0.021 ± 0.016 


0.013 


0.020 ± 0.024 


0.008 


0.015 ± 0.018 


0.039 


0.060 ± 0.076 


0.028 


0.074 ± 0.108 


0.037 


0.044 ± 0.041 


0.026 


0.037 ± 0.034 


0.013 


0.017 ± 0.014 


0.007 


0.011 ± 0.014 



Table 8. Resultant KS probabilities* concerning emission-line luminosities normalized by i/L u (25 /xm). 

Line Sl to tai vs. SI. 5s Sl to tai vs. S2 to tai SI. 5s vs. S2 to tai 

L([0 ni]A5007)/z/Z, [ ,(25/mi) 3.513X10" 1 1.800X10" 1 1. 280x10^ 

L([0 iii]A4363)/i/L„(25/zm) 2.916X10" 1 2.474xl0~ 2 3.038X10" 1 

L([Nc m]A3869)/i/X l/ (25/xm) 4.313X10" 1 2.450x10-* 6. 721x10^ 

L([Fe vii]A6087)/t/£„(25/mi) 4.983xl0~ 2 1.801xl0~ 3 4.809 xlO" 2 

These values means the probabilities of two sample drawn from the same parent population. 
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log z 



Fig. 1. — Frequency distributions of the redshift for each class of Scyfert nuclei. The number of objects for each types is written in 
parenthesis. Note that the data of a nearby SI. 5, NGC 3031, is not included in this figure because its redshift is negative. 




40 42 44 

\ogv L v (25jum) (ergs s _1 ) 

Fig. 2. — Same as figure 1, but for the IRAS 25 fim luminosity. 
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Fig. 3. — Same as figure 1, but for the IRAS 60 fim luminosity. 
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Fig. 4. — Compiled objects are shown in the diagnostic diagrams proposed by Veilleux and Osterbrock (1987). The diamonds are NLSls, 
the squares are BLSls, the circles are SI. 5s, the triangles pointing to the left are S2r,blrs, those pointing to the right are S2hblrs, and those 
pointing to the upper side are S2~. Red color means Sis, blue color means SI. 5s, and green color means S2s. For comparison, compiled H 
II systems are also plotted in this figure. The plus signs are blue compact galaxies (Izotov et al. 1994), the crosses and the "Y" signs are 
extragalactic H II galaxies (McCall et al. 1985; van Zee et al. 1998), and the asterisks are H II galaxies (Masegosa et al. 1994). (a) Diagram 
of [O in]A5007/H/3 versus [N n]A6583/Ho. (b) Diagram of [O ni]A5007/H/3 versus [S n]AA6717,6731/Ho. (c) Diagram of [O iii]A5007/H/3 
versus [O l]A6300/Ha. 
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Fig. 5. — Diagram of H7/H/3 versus Ha/H/3. The symbols are the same as those in figure 4. Theoretically predicted emission-line ratios 
assuming the case B and dust extinction of mag < Ay < 5 mag are also plotted by a solid line. 




log ([O l]A,6300 / [O lll]A.5007) 



Fig. 6a. — Frequency distributions of the various emission-line flux ratios, (a) Flux ratio of [O l]A6300/[O in]A5007. (b) Flux ratio of [O 
n]A3727/[0 m]A5007. (c) Flux ratio of [O i]A6300/[O n]A3727. (d) Flux ratio of [O m]A4363/[0 m]A5007. (e) Flux ratio of [S n]A6717/[S 
Il]A6731. (f) Flux ratio of [O i]A6300/[S n]AA6717,6731. (g) Flux ratio of [O n]A3727/[S n]AA6717,6731. (h) Flux ratio of [S n]AA6717,6731/[0 
m]A5007. (i) Flux ratio of [O i]A6300/[N o]A6583. (j) Flux ratio of [O o]A3727/[N o]A6583. (k) Flux ratio of [N n]A6583/[0 ni]A5007. (1) 
Flux ratio of [S n]AA6717,6731/[N n]A6583. (m) Flux ratio of [Ne m]A3869/[0 m]A5007. (n) Flux ratio of [Ne m]A3869/[0 n]A3727. (o) Flux 
ratio of [Ne v]A3426/[0 n]A3727. (p) Flux ratio of [Fe vn]A6087/[O m]A5007. 
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log([OII]A3727/[OIII]A5007) 



Fig. 6b. — continued. 




log([OI]A,6300/[OII]A,3727) 



Fig. 6c. — continued. 
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log ([O IIIJA.4363 / [O IIIJA5007) 



Fig. 6d. — continued. 
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Fig. 6e. — continued. 
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Fig. 6f. — continued. 
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Fig. 6g. — continued. 
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Fig. 6h. — continued. 
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Fig. 6i. — continued. 
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log ([O II] A.3727 /[Nil] Aj6583) 



Fig. 6j. — continued. 
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Fig. 6k. — continued. 
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Fig. 6l. — continued. 
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Fig. 6m. — continued. 
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Fig. 6n. — continued. 
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Fig. 60. — continued. 
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Fig. 6p. — continued. 
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Fig. 7. — Frequency distributions of the various emission-line flux ratios. Since the number of objects in each sample is small, the histograms 
are shown only for the Sl to talS, the SI. 5s, and the S2 total . (a) Flux ratio of [O n]A7325/[0 n]A3727. (b) Flux ratio of [S ni]A9069/[S 
ll]AA6717,6731. (c) Flux ratio of [N l]A5199/[N ll]A6583. (d) Flux ratio of [N ll]A5755/[N ll]A6583. (e) Flux ratio of [Ar in]A7136/[0 m]A5007. 
(f) Flux ratio of [Fe x]A6374/[0 m] A5007. (g) Flux ratio of [Fe xi]A7892/[0 m] A5007. 
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Fig. 8a. — Various emission-line flux ratios are plotted as functions of redshift, IRAS 25 fim luminosity, and IRAS 60 (im luminos- 
ity. The symbols are the same as those in figure 4. (a) Flux ratios of [O l]A6300/[O m]A5007, [O ll]A3727/[0 m]A5007, [O l]A6300/[O 
n]A3727, [O in]A4363/[0 m]A5007 (=floin), [S n]A6717/[S n]A6731 (=flsn), and [O i]A6300/[S n]AA6717,6731. (b) Flux ratios of [O 
n]A3727/[S n]AA6717,6731 [S n]AA6717,6731/[0 m]A5007, [O i]A6300/[N n]A6583, [O n]A3727/[N n]A6583, [N n]A6583/[0 m]A5007, and [S 
n]AA6717,6731/[N o]A6583. (c) Flux ratios of [Ne m]A3869/[0 m]A5007, [Ne m]A3869/[0 n]A3727, [Ne v]A3426/[0 n]A3727, [Fe vo]A6087/[O 
in]A5007, [O n]A7325/[0 n]A3727 (=R n), and [S m]A9069/[S n]AA6717,6731. (d) Flux ratios of [N i]A5199/[N n]A6583, [N ll]A5755/[N 
n]A6583 (=R N n), [Ar m]A7136/[0 m]A5007, [Fe x]A6374/[0 m]A5007, and [Fe xi]A7892/[0 m]A5007. 
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Fig. 8b. — continued. 
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Fig. 8c. — continued. 
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Fig. 8d. — continued. 
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Fig. 9. — Frequency distributions of the luminosities of the [O m]A5007 emission, the [O m]A4363 emission, the [Ne m]A3869 emission and 
the [Fe vn]A6087 emission for the Sltotal s i the SI. 5s, and the S2 tota is. These luminosities are normalized by the IRAS 25 (im luminosity. 
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FlG. 10. — Diagram of L([0 m]A5007) versus L([0 ll]A3727). The symbols are the same as those in figure 4. 



o 
o 



O 10' 



CO 
CD 



CD 



w 



10" 





' 






■ 


, , | 


' 


. 


A 














" 




A 


<1 










" 




<1 

< 
< 

O 




O 










~~ 


< A A ^ D 






<1 


O 




— 




A 






o o 














4 




A □ 




. 




A 

A^ 


A 


A 




< < 




. 




C 




o 


> „ D 












9 


A 


A OO 








o 

A 


A 


A 


m 


<>o 








Cfl 


A 


A °A @ □ 




■ 




< 




o 


n 


* ,A A 








<i 


□ 




Ai 


> A ° ° □ " 

a < " o & i : d> 
o o° o o 


O 

D 






<l 
















< 




D 




'■ 












r-, O 


















<o 


■ 












. . 1 




■ 



10" 



10 L 



vF v {25jum) / vFJieOjum) 



Fig. 11. — Diagram of the emission-line flux ratio of [S ll]AA6717,6731/[0 m]A5007 versus that of vF„(25 /tm)/i/F„(60 fim). The symbols 
are the same as those in figure 4. 
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Fig. 12.— Calculated emission-line flux ratios of [O m]A4363/[0 m]A5007 (Rom), [Ne in]A3869/[0 m]A5007 and [Fe vn]A6087/[O m]A5007 
are shown as a function of the gas density. Here, the Lyman optical depth is assumed to be rgi2 = 0.1. The dotted lines, the dashed lines, 
and the solid lines denote the models adopting U = 10 — 3,s , 10 — 25 , and 10~ 1,5 , respectively. The blue lines, the green lines, and the red lines 
denote the models adopting the mctallicity of half the solar, the solar, and twice the solar one, respectively. 
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Fig. 13. — Calculated emission-line flux ratios of [Ne m]A3869/[0 m]A5007 and [Fe vn]A6087/[O m]A5007 shown as a function of the 
ionization parameter. In these models, njj = 10 60 cm -3 is assumed. The solid lines denote the models with the metallicity of half the solar 
value, the dashed lines denote the models with the solar metallicity, and the dotted lines denote the models with the metallicities of twice the 
solar value. 
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Fig. 14. — Same as figure 12, but for rgi2 = 0.01. Only the models with solar metallicity arc shown. 
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Fig. 15. — Same as figure 12, but for rgi2 = 1.0. Only the models with solar metallicity arc shown. 



